William I. Fine Theoretical Physics Institute 
University of Minnesota 



FTPI-MINN-08/36 
UMN-TH-2718/08 

OO ' 

September 2008 

' 

o 

Q 

a> ! Diresonance in production and scattering of heavy 

^ mesons 
PK 

S. Dubynskiy 

1 / School of Physics and Astronomy, University of Minnesota, 
^ ! Minneapolis, MN 55455, USA, 
^ ! M.B. Voloshin 

William I. Fine Theoretical Physics Institute, University of Minnesota, 
cn ; Minneapolis, MN 55455, USA 

g i and 

Institute of Theoretical and Experimental Physics, Moscow, 117218, Russia 



X 



Abstract 

We consider the production and scattering amplitudes of heavy mesons in a situa- 
tion, where there are two closely spaced narrow resonances, which structure we refer to 
as a diresonance. Assuming strong overlapping of the resonances coupled to common 
channels, it is found, using the unitarity and analyticity constraints, that the produc- 
tion amplitudes by a weak source should have similar behavior with energy in different 
channels. In particular the ratio of the coefficients for each pole contribution to the 
production amplitude is fixed at —1. 



The spectroscopy of resonances near the charm threshold attracts a considerable renewed 
interest. Recent experimental data not only provide evidence of new states but also suggest 
that the nature of well known resonances merits a reexamination at the level of fine details. 
The new exotic charmonium-like states [T] considerably expand the standard spectrum of 
charmonium, and challenge us for a better understanding of the strong dynamics. Further- 
more the apparent significance of multiquark states at the onset of open charm threshold may 
also impact the properties of the known resonances and may be instrumental in resolving 
some long-standing puzzles. One such puzzle is related to an inconsistency between BES and 
CLEO results for the production and decays of ?/>(3770) resonance, in particular the fraction 
of its decays into non-DD states[2J[3]. Recently BES Collaboration has reanalyzed their data 
on e + e~ annihilation in the energy region between 3.700 and 3.872 GeV They reported 
observation of an anomalous line-shape behavior of the cross section which is inconsistent 
with the presence of only one simple z/>(3770) resonance in this energy region. It is claimed 
that this anomalous behavior could be better understood in terms of two resonances near 
the cm. energies of 3.764 GeV and 3.779 GeV. This result could violate the conventional 
interpretation of ?/>(3770) as being a dominantly l 3 Di charmonium state with an admixture 
of 2 3 Si, and clearly suggests a more complicated structure, possibly including strong dy- 
namics of the D meson pairs near the threshold. Namely, a diresonance structure may arise 
from existence of both a charmonium state and a 'molecular' DD threshold resonance. 

Apriori one would expect that the nature of each of the two individual peaks could be 
studied by further exploring their relative coupling to various channels. However the pur- 
pose of the present paper is to argue that this standard method, applicable to sufficiently 
widely separated resonances, is unlikely to be applicable for a strongly overlapping pair of 
resonances, such as the one indicated by the BES data, i.e. when the splitting between 
the positions of the resonances is comparable with their widths, and all of these parameters 
are small in a typical energy scale for the process. (It is natural to call such a structure 
as diresonance.) Using unitarity and analyticity constraints, we find, under the simplest 
assumptions, similar to those involved in the standard Breit-Wigner treatment of a single 
resonance, that the two states are necessarily strongly mixed, and that various final states 
produced e.g. in the e + e~ annihilation have the behavior of the production cross section in 
the diresonance region proportional to one another. Furthermore, if a diresonant production 
amplitude is written as a linear combination of two poles, then in the limit where the dires- 
onance parameters, the two widths and the splitting between the poles, can be considered 
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as small, the relative factor between the two poles is necessarily equal to —1. 

In practical terms this universal behavior implies that it would be problematic to dis- 
entangle experimentally the underlying origin of the states in a diresonance complex. In 
particular, even if the suggested by the experiment structure around 3.77 GeV, originates 
from an overlap of charmonium and molecular states, the mixing between them effectively 
erases any difference in their experimental signature. 

We start our discussion of a diresonance structure with a simple case of just one scat- 
tering channel. The radial part of the wave function of the particle with mass m with 
arbitrary complex energy E and orbital momentum I = (here we consider S'-wave motion 
for simplicity) has the following form at large distances r: 

R= - [B*(E)e ikr + B(E)e~ ikr ] , k = V2mE , (1) 

where the coefficients of the incoming and outgoing waves are related by complex conjugation 
due to the requirement that the wave function is real at real negative E. In the familiar 
case of a single resonance at E r = Eq — iT/2, with Eq and T being the position and the 
width of the resonance, one has B(E r ) = 0, which ensures that the wave function of the 
resonant state vanishes at spatial infinity [5]. One can then expand the function B(E) near 
the position of the resonant level E r as 

B(E) = (E-E + l -T)b, (2) 

with b being a smooth function of energy, i.e. b changes on a scale much larger than the 
resonance width T. The scattering S- matrix element is then found as follows 

where 5^ is a nonresonant phase which is a smooth function of energy defined as exp(2i5^) = 
b*/b. 

This standard Breit-Wigner treatment of a single resonance can be readily extended to 
the case of a diresonance, i.e. in the situation when the scattering amplitude has two closely 
separated poles at E\ —iT\/2 and E 2 — iT 2 /2 with both widths Fx and T 2 and the difference 
Ei — E 2 being considered as 'small'. The expansion of the coefficient B(E) having two zeros 
in the diresonance region is obviously given by 

B(E) = (Ai+t 1 i)(A 2 + t l2 )b, (4) 
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where the notation is introduced A a = E — E a , 7a = T a /2 (a = 1, 2), and, similarly to Eq.p]), 
the coefficient b is a slowly varying function of the energy, which can be approximated by 
a constant on the energy scale of the diresonance region. The corresponding expression for 
the S matrix element then takes the form 

_ (Ai-*7i)(A 2 -?7 2 ) , (o), f5l 
5 -(A 1 + * 7l )(A 2 + * 72 ) eXp(2 ^ )' & 

Let us consider now the amplitude A(E) for production of the scattering state in the 
diresonance region by a point-like source. The production process is assumed to be weak, so 
that it is sufficient to consider only the lowest order in the coupling to the source. The energy 
dependence of such amplitude is proportional to the inverse of the coefficient of the incoming 
wave in the wave function ([1]): A(E) = g/B(E) with g being a real (at real E) smooth 
function of energy. Indeed, according to the familiar u if)(0) rule" the absolute value of the 
amplitude is proportional to ip(0), provided that the wave function is normalized to a fixed 
amplitude at infinity, R = (1/r) sin(/cr + 8), which implies the relation cx 1/\B(E)\. 

On the other hand, according to the Watson's theorem, the phase of A is given by 8, i.e. the 
phase is that of 1/B(E). Using this relation we readily find an analytical formula for the 
production amplitude in the diresonance energy region 

ME) = , A . Wfl • (6) 

v ; (A 1 + z 7l )(A 2 + « 72 ) v ; 

The latter expression for the diresonance production amplitude when written as a sum 
over two resonances: 

A{E) W) , (i.i) (7) 

E 1 - E 2 + ry 2 - r/i \ A 1 + 1 7l A 2 + 1 72 J 

tells us that the relative phase between the two resonance factors has to be equal to it and 
the coefficients of the pole factors should be the same. 

As a simple cross check we considered a toy model with the scattering of a particle with 
mass m in a central potential with two Gaussian barriers: 



V ( r ) = o 2 X 1 kl eX P 



(x — X\) 2 ~\ i , [ (x — x 2 ^ 2 



Wi 



h 2 exp 



w 2 



(8) 



where x = r/ro is a dimensionless ratio of the distance r to an arbitrary scale ro and 
the parameters hi, also dimensionless. The distance scale r also sets the scale 

Eq = (2mrQ)~ 1 for the energy of the particle. We calculated numerically the dimensionless 
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"production amplitude" as the value of the wave function at the origin, ^(0), at energy E, 
provided that at large r the wave function is normalized to a wave with unit amplitude. 
We found that every time the parameters of the potential hi, tuned in such a way 

that a diresonance structure appears, the production amplitude is closely approximated by 
the expression (j7j). We illustrate this behavior in Fig. 1 for one specific set of parameters 
(hi = 5.5, hi = 0.94, xi = 1.75, X2 = 6.92, w\ = 0.86,^2 = 0.45). The fit curve shown in the 
plot corresponds to a constant ratio g/b in Eq.fJTJ, so that the relative factor between the 
two poles exactly equals —1. If the fit is relaxed and this relative factor is also treated as a 
fit parameter, we find the best approximation for it as —0.90 + 0.02i. 
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Figure 1: The dimensionless "production cross section" |A(i?)| 2 vs the energy E (in units of 
Eq) in the diresonance region in a toy model with potential scattering. The circles are the 
numerical data and the curve is the fit by the formula in Eq.(J7J). 



It may appear that the rigid constraint on the relative contribution of the two single- 
resonance factors in a diresonance complex is a limitation of the considered situation with 
one scattering channel. For this reason we proceed to discussing scattering in two coupled 
channels, where we find that the same constraint on the relative strength of the two pole 
factors applies in each of the channels, so that in fact the production of the final states in 
the two channels is described by the same energy dependence with the only free factor being 
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the ratio of the overall yield of each of the final channels. 

In this discussion we consider two channels a and b and consider the S matrix for the 
scattering and the production by a localized source: 

Saa i 

S = S ba S bb %A h , (9) 
\iA a %A h 1 / 

where A a ^ are the production amplitudes for each of the channels by a weak localized source, 
so that the quadratic in the strength effect of the source in the diagonal (33) element of the 
S matrix can be neglected. In the zeroth order in the production amplitudes we find that 
the general solution to the unitarity conditions for the two-channel scattering matrix with a 
diresonance singularity has the form 

2i (Ai 72 + A 2 7i) rj m r) n ~ 



Q 



5mn (A 1+ * 7 i)(A 2 + *7 2 ) J 



exp 



;<(*£> + #>)] ; m,n = a,b, (10) 



where S a °l is the nonresonant scattering phase in the corresponding channel and the real 
factors rj a and r\ satisfy the condition rj^ + rjl — 1. It can be noted that in the single 
resonance case these factors are determined by the corresponding branching fractions for 
the resonance: rj 2 = T n /T. In the diresonance case we do not find a simple direct relation 
of these factors to the individual width parameters I\ and T 2 . However these factors can 
still be interpreted in terms of the branching ratios (for the diresonance complex) in the 
sense that rfc gives the probability of the branching of the scattering into the corresponding 
channel. It is also quite clear that setting rj to zero in one channel and rf = 1 in the other, 
returns us to the previously discussed case of a single channel. 

The unitarity relation for the matrix (jH]) in the first order in the production amplitudes 
A atb then yields these amplitudes in the form 

A n = , A —. ^w" - . — r exp(^ 0) ); n = a, b, (11) 
(Ai + z7i)(A 2 + «7 2 ) 

where the real smooth factor fi characterizes the strength of the coupling of the source. 

The expression ( TTTI) implies that the single- channel behavior of a diresonance production 
amplitude also holds for multiple channels. Namely, when expressed as a linear combination 
of two poles the relative coefficient between the two pole terms is necessarily equal to — 1 
as shown in Eq.^7]), and the relative yield in each channel is determined by the branching 
factor rj 2 . 
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One can readily notice that the relative factor —1 is a trivial consequence of the produc- 
tion amplitude behaving as (AE)~ 2 away from the diresonance region. Clearly, in order to 
invalidate such a behavior, one would have to introduce in the dependence of the coefficient 
g/b an energy scale comparable to the parameters (Ei — E 2 ), I\ and T 2 . Thus our conclusions 
are applicable in the situation, which we refer to as a diresonance, where these parameters 
are small in the typical scale in the problem. In this respect the assumption is quite similar 
to the familiar Breit-Wigner approximation for a single resonance, where the resonance has 
to be considered as narrow, i.e. with a small width. In the diresonance case it is also the 
splitting between the poles, which has to be "narrow" in addition to the width parameters, 
for our approximation to be valid. 

We believe that our consideration of a diresonance structure is quite generic and may 
be applicable to the suggested by experiment [2] structure in the e + e~ cross section near 
3.77 GeV, and possibly to other similar structures. A specific detailed application of the 
discussed diresonance properties to the production of D meson pairs in the ^(3770) region 
should also include the P-wave kinematics with different thresholds for the pairs of neutral 
and charged mesons, as well as the Coulomb effects. Such analysis can be done along the 
lines presented in Ref. [6], if more detailed data become available. At this point we can only 
remark that the two-pole fit to the data [I] does not contradict the expression (JTj) for the 
production amplitude. However the error range is still too large for any further conclusions 
to be drawn. 

This work is supported in part by the DOE grant DE-FG02-94ER40823. The work of 
S.D. is supported in part by the Stanwood Johnston grant from the Graduate School of the 
University of Minnesota. 
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